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Abstract 

In this paper, the distribution of mass and kinetic energy in the cold region of the thermal 
neutron induced fission of ^^^U, 235-y ^^^^^ 239p^^ respectively, is interpreted in terms of nucleon 
pair-breaking and the Coulomb interaction energy between complementary fragments (Coulomb 
effect). The fission process ends in the scission configuration, in which complementary fragments 
1 and 2, having proton numbers Z\ and Z2, and mass number A\ and A2, only interact each other 
with electrostatic repulsion, by what each one acquires kinetic energy E\ and £"2, respectively. 
However, before reach detectors, they emit n\ and n^ neutrons. Thus, the final detected fragments 
have mass numbers m\ (= A\ — n\) and m^ (= Ai —n2), and kinetic energies ei (= Ei[l — ni/Ai]) 
and 62 (= -E'2[l — "'^2/^2])! respectively. In order to avoid the erosive consequences of neutron 
emission, one studies the cold fission regions, corresponding to total kinetic energy (TE) close to 
the maximum available energy of the reaction (Q). Contrary to expected, in cold fission is not 
observed high odd-even effect in mass number distribution. Nevertheless, the measured values are 
compatible with higher odd-even effects on proton or neutron number distribution, respectively. In 
addition, in cold fission, the minimal total excitation energy (X) is correlated with the "Coulomb 
energy excess" (SC), which is defined as the difference between C (the electrostatic interaction 
energy between complementary fragments taken as spherical in scission configuration) and Q. 
These Coulomb effects increase with the asymmetry of the charge fragmentations. In sum, the 
experimental data on cold fission suggest that scission configurations explore all the possibilities 
permitted by the available energy for fission. 
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I. INTRODUCTION 

Among the most studied properties of the nuclear fission of actinides are the distributions 
of mass, charge and kinetic energy associated to complementary fragments. However, R. 
Brissot et al. [l| observed that the neutron emission from the fragments, occurred before 
they reach the detector, generates a structure on the final kinetic energy distribution as 
a function of mass, non-existent on primary fragment distribution. Naturally, neutron 
emission also erodes odd-even effects on neutron number distribution and consequently on 
mass number distribution associated to fragments. For these reasons, cold fission regions, 
corresponding to high kinetic energy associated to fragments that do not emit neutrons, 
are studied. In that sense, in this paper, mass, charge and kinetic energy distribution in 
the cold regions from thermal neutron induced fission of ^'^^U, ^^^U and ^'^^Pu, namely 
^^^U(nt/i,f), ^^^U(nt/i,f) and ^^^Pu(nt/i,f), respectively, is interpreted in terms of nucleon pair 
breaking and Coulomb interaction energy between fragments in scission configuration. 

Let us suppose that the complementary fragments have kinetic energies Ei and E2, 
referred to primary mass numbers Ai and A2, and proton numbers Zi and Z2, respectively. 
After the neutron emission, complementary fragments 1 and 2 end with mass numbers 

nil = Al — ni 

and 

1712 = A2- n2 

where rii and n2 are the numbers of neutrons emitted by the fragments 1 and 2, respectively. 
The corresponding values of the final kinetic energy associated to those fragments are 



and 



e.^B.|l-^|. 



In the case that only one of the complementary fragments is detected, which is what happens 
in the Lohengrin spectrometer at the High Flux Reactor (HFR) of the Laue-Langevin 
Institute (ILL) in Grenoble [l|, subscripts are omitted. For a given primary mass number 
{A), the kinetic energy (E) distribution is characterized by the corresponding average, 



E{A), and standard deviation, aE^A). 

Using the Lohengrin spectrometer in the HFR of ILL in Grenoble, R. Brissot et. al. 
measured e(m) y o"e(m) values in ^^^U(nt/i,f). Their experimental curve ae{m) showed a 
peak around m = 110, as it had been predicted by a Monte Carlo simulation performed 
before that the experiment was carried out. In this simulation it was demonstrated that 
that broadening e{m) distribution is due to the neutron emission and the sharp fall on 
E[A) in the mass region around m = 110. For each mass number m,the kinetic energy 
distribution is the result of superposition of distributions associated to fragments having 
A>m\^. 

In the final mass region 106 < m < 112, e falls sharply, about 4 MeV per mass 
number. Let us suppose that the E behaviour is similar to that of e, and (7e{A) is flat 
(equal to 5 MeV, for example). Therefore, the average kinetic energies associated to 
fragments with A = 112 is 8 MeV lower than that of fragments associated to A = 110. 
Fragments with A = 110, that do not emit neutrons, have the high values of their 
distribution of primary energies (E); and fragments associated to A = 112, that do emit 2 
neutrons, have the low E values of their corresponding primary distribution, whose average 
is already 8 MeV lower than that associated to A = 110. Both fragments, one with primary 
mass numbers 110 y and the other with 112, respectively, having E values with average 
differences of around 8 MeV, end with m = 110. The superposition of these dispersed 
values contributes with the widening of ae{m) = 110. 

D. Belhafaf et. al, repeating the experiment on ^^^U(nf/j,f) made by R. Brissot et. 
al. l|, observed a broadening of the distribution of e around m = 126 J2|. Moreover, for 
^^^U(ntft,f), those authors observe a broadening of the distribution of e around m = 124. 
Their simulation did not reproduce these experimental results 



a 



For ^^^U(nj/i,f) H. Faust et. al. [3|, using theoretical models on fission dynamics, 
calculated the distribution of i? as a function of A, but not around A = 124. However 



the D. Belhafaf et. al. experimental results for both, ^^'^U(nt/i,f) and ^^^U 
reproduced by a Monte Carlo simulation performed by de M. Montoya et. al. 



ntfe,f), were 
4, l5| . These 



authors shows that the broadening on the distribution of e is generated by the superposition 
of the distributions of E corresponding to A > m. For the case of ^^^U(nt/i,f), e falls from 
88 MeV referred to m = 125, to 84 MeV referred to m = 126, while in this region the 



a 



tendency of e is to increase with m [2|. See Fig. [T] It is reasonable to assume that for 
A near but higher than 126, E has a similar behavior. Thus, the superposition E values, 
corresponding to primary mass numbers A > 126, having dispersive behavior, contribute to 
the broadening of e{m) distribution around m = 126. 

II. NEUTRON EMISSION AND ODD-EVEN EFFECTS IN COLD FISSION 

The odd-even effect on the fragment mass number distribution {5A) is defined by the 
relation: 

where Y^e and Y^^ are yields associated to fragments having even and odd mass numbers, 
respectively. 

Similarly, we define the odd-even effects 6N and 6Z on the distribution of neutron 
number (A^) and proton number [Z), respectively. 



For low energy fission of actinides, the odd-even effect on the proton number distri- 
bution of the fragments, previously found to A. L. Wahl et. al. |6[ and G. Mariolopoulos et. 
at. J7|, suggested that must exist high odd-even effects on the mass number distribution, 
but they are eroded by neutron emission. 

To study the mass number distributions undisturbed by the emission of neutrons, C. 
Signarbieux et. al. chosen high kinetic energy windows associated to fragments from 
^^^U(nt/i,f), ^^^U(nt/i,f) and ^^^Pu(nt/i,f), respectively. The number of events analyzed were 
1.5x10^, 3x10^ and 3.2x10^, respectively. These authors used the difference of time of 
flight technique to separate fragment mass numbers, with solid detectors to measure the 
fragment kinetic energy. The experiment was conducted in the HFR of ILL in Grenoble |8|,|9|. 

Although 5A is not zero, it is not as high as was expected. However, one may show 



that this small SA is compatible with nucleoli pair breaking and higher odd-even effects in 
the distribution of proton or neutron numbers. 



H.-G. Clerc et. al. 10|, for E = 114.1 MeV, obtained the mass number distribution 



presented in Fig. ^ from which one calculates 6 A = 30.14 %. The mass number distribu- 



tion in the window of 
is presented in Fig. [3] 



ight fragment kinetic energy 113.5 < E < 114.5 MeV, in ^^^U(ni/i,f), 
9|. From this distribution one gets 6 A = 32.18 %. 



In order to interpret values of 6A, 6Z and 6N, lets first define, for each A, Zq as 
the light fragment proton number that corresponds to maximal available energy of the 
reaction (Q) Suppose that nucleon pair-breaking occurs, but fragment couple reaching 
maximum total kinetic energy still corresponding to Zq; therefore, for ^'^^U(n(/t,f), 6Z = 66.6 
% and 5N = 33.3 %. This result is coherent with the stereographic projection of the values 
obtained by U. Quade et. al. ll| for lower energy windows, as it is presented in Figs. |4]and[5l 



For total kinetic energy TE > 204 MeV, only ^°^Mo {Z = 42, N = 62) survives 
therefore 6 A = 1. 



The mass number distribution for the kinetic energy window 113.5 < E < 114.5 
MeV, for light fragments from the reaction '^^^U{nth,i), is presented in Fig. [6l from which 
one calculates 6A = 15.3 %. For E = 113.7 MeV, W. Lang et. al. obtained the distribution 
presented in Fig. [3, from which one calculates SA = 23.82 %. For E = 108 MeV, W. 



Lang et. al. obtain SA = 0, SZ = 35 % and SN = 8 % 13] • Suppose that proton number 



correspond to Zq, then one calculates SZ = 66.6 % and SN = 33.3 %. This result is 



coherent with the stereographic projection of the values obtained by W. Lang et. al. 12] 
for windows of lower energies, as it is shown in Figs. [8] and [91 



For total kinetic energy TE > 203 MeV, the surviving mass numbers are A = 104- 
106, among which there are 16 even and 8 odd mass numbers; therefore SA = 30 %. 
For TE > 200 MeV J. Trochon et. al. [Ij] measured mass number and proton number 
distributions. From their results one obtains SA ^ 20 %; SZ ^ 60 %, and SN ^ 60 %. 



The mass number distribution in the window of hght fragment kinetic energy E 



112 MeV, in ^^^Pu(n(/i,f), obtained by C Schmitt et. ai, is presented in Fig. [TO] 15|. From 



this data one calculates SA = 3.7 %. For this energy they obtain SZ = 15 %, and SN = 10 



% 15|. For E > 119 MeV, one obtains the distribution presented in Fig. [TTl from which 
one calculates 6 A = 7.6 %. Suppose that for E' > 119 MeV, the proton numbers correspond 
to Zq, therefore 6Z = 66.6 % and 6N = 33.3 %, which are coherent with the stereographic 
projection of the values obtained by C Schmitt et. ai, for lower energy windows. See Figs, 
and I 



For TE > 210 MeV, one obtains masses 104 < A < 107 , with the higher yield cor- 
responding to A = 106. 

For lower kinetic energy windows, one obtains SA = 0. This result does not imply 

that 6Z y 5N values are zero. M. Montoya et. al. using a combinatorial analysis 
demonstrated [9| that if there is one and only one nucleon (proton o neutron) pair breaking, 
then 

l + 6A = 6N + 5Z. (1) 



H. Nifenecker et. al. 16|] developed other combinatorial analysis that produces the same 
relation. 

III. SHELL AND COULOMB EFFECTS ON COLD FISSION 



In cold fission from ^^'^U(nj/j,f), ^^^U(ni/i,f) and ^^^Pu(nt/i,f), shell and Coulomb effects 
on mass number and kinetic distribution were observed 



'u(n 



To describe "Coulomb effect" let's first define the following curves depending of A: 
i) Q, the smoothed curve referred to maximal value of available energy of the reaction (Q), 
as a linear curve in the mass region from the lightest A to that at which ends the linear 
trend and, for higher mass numbers, a curve following the approximately horizontal trend 
of Q; ii) C, the electrostatic interaction energy between complementary spherical fragments 
corresponding to Zo,with surfaces separated by 2 fm (the result is a step function of mass 
number); and iii) K., the maximal value of total kinetic energy, as the lowest of the ten 



highest values of total kinetic energy, defined as presented by M. Montoya for ^^^U(nift,f), 
235u(ni,„f) and ^^^Fn{nth,i), respectively pj. 

Now let's define the smooth excess of electrostatic energy, as a function of A, by the 
relation 

6C = C-Q (2) 

and the smoothed values of minimal excitation energy, as a function of A, by the relation 

X = Q-K (3) 

The curves of C, K, and Q values for ^'^^U(ni/i,f) are presented in Fig. UM 

The fiuctuations of K in ^^^U(nt/i,f), ^^^U(ni/i,f) and ^^^Pu(nt/i,f) are interpreted as 



effects of the step variation of the electrostatic energy as a function of Zq 17|, [iSj 



It is found that the smoothed values 6C and X are correlated. This result suggests 
that, at scission, a higher excess of electrostatic interaction energy produces a higher 
deformation of fragments, therefore a higher X and A then lower K- values. 

For the mass regions 78 < A < 100 in ^^^V{nth,i); 80 < A < 102 in '^^^V{nth,i); 
and 90 < A < 108 in ^^^Pu{iath,i)'- SC and X are decreasing functions of A. See Figs. [T5l 
[I6]and[171 

In ^^^U(ni/i,f), the higher Q- values, are around 204 MeV, which occurs in the mass 
region 100 < A < 106, which corresponds to the lowest values of 6C. See Fig. [151 For 
complementary fragments {Z = 50, N = 80) and {Z = 42, A^ = 62), the K value reach the 
respective Q value. 

As it is presented in Fig. [111 in the mass region 100 < A < 106, in ^^^U(n(/j,f), 6C 
well as X are very approximately zero. 

In^^^U(nift,f), the lowest X value occurs in the mass region 100 < A < 110. See 
Fig. [13 

8 



For ^^^U(n(/i,f), '^^^l]{nth,i) and ^^^Pu(n(/i,f), the value of SC is not zero. At least 
one of the fragments must be deformed and spends energy for that. However, molybdenum 
{Z = 42) with A^ = 60, 62 y 64; and zirconium {Z = 40) with N = 60, 62 y 64; have prolate 
geometries in their ground states, in addition they are transitional and soft [9]. At scission, 
these fragments are deformed without excitation energy, and the electrostatic energy is 
equal to the Q- value. 



For -U(„,a). -U(n,a)„and -P„(n,a). the fluctuations of ,„a^ number yeld 



are anti-correlated to 6C p, l7| 



For ^^^U(n(/i,f), in the light fragment kinetic energy window 113.5 < E < 114.5 MeV, the 
mass yield peaks toward a maximum around A = 90 and 100, respectively, which correspond 
to minimal values of 5C |9|, [l6|. For ^^^U(ni/i,f), in the light fragment kinetic energy window 
113. 5< E <114.5 MeV, the mass yield peaks toward a maximum around A = 90 y 102, 
respectively, which correspond to minimal values of 5C. For ^'^^Pu(n(;j,f), in the light 
fragment kinetic energy region i? > 118 MeV, the mass yield peaks toward a maximum 
around A = 84, 89, 96, 100 and 106, respectively, which correspond to minimal values of 6C. 

For ^^^U(ni/i,f) as well as for ^^^U(n(ft,,f), fori? > 106, SC increases sharply with A 
and, therefore, K decreases sharply too. For the reaction ^^^U(nj/t,f), 5C decrease for 
A> 110. 

The fluctuations of X and 6C are correlated and they increase with asymmetry of 
fragmentations. This is because the step-variation of the electrostatic interaction energy 
between fragments produced, by changing a unit of proton number, is more pronounced in 
asymmetric region, as one may in Fig. [171 corresponding to ^'^^U(nt/i,f). 

The electrostatic interaction energy between complementary fragments constitutes a 
barrier to fission in the channel of the corresponding charge and mass. The probability 
to produce a couple of fragments will be higher for lower difference between this barrier 
and the available energy. Thus, for a same Q-value, between neighbor mass numbers, the 



fragments with higher asymmetry of charge will reach higher values of kinetic energy. This 
is equivalent to say for TE close to the maximal value the prevalent fragmentations will be 
those with higher asymmetry. The Coulomb effects in cold fission will be present also in 
proton number yields corresponding to odd proton number. The higher yield corresponds 



to higher asymmetry, as it was shown by W. Schwab et. al. in the case of ^^'^U(nt/i,f) 19|. 



IV. CONCLUSIONS 

Neutron emission has erosive consequences on fragment mass and kinetic energy distri- 
bution in ^^^U(nt/i,f), ^'^^U(nj/j,f) and ^^^Pu(nj/j,f). To avoid those disturbances, the cold 
fission, i.e. very low excitation region, was studied. Contrary to expected, in cold fission, 
nucleon pair breaking still exists. Nevertheless at the maximal total kinetic energy only 
even-even fragments survive. For lower energies, zero odd-even effect on mass distribution 
does not exclude non zero odd-even effects on proton or neutron distributions. 

In addition, in cold fission. Coulomb effects on charge, mass and kinetic energy dis- 
tribution are present: for two fragmentations with equal Q-value, the higher yield 
corresponds to the higher charge asymmetry. Moreover, Coulomb effects are presented as 
fluctuations of the maximal kinetic energy value as a function of mass. This effect is due to 
the fact that, each five units of mass number, there is a change of the proton number that 
maximizes the Q-value, which produces a step variation of electrostatic interaction energy 
between complementary fragments with the same shapes. A similar effect on mass number 
yield curve is observed. Due to the characteristics of electrostatic interaction, those effects 
are more notorious in higher asymmetry region. Finally, results from several experiments 
in cold fission suggest that scission explore all configurations permitted by the Q-value. 
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FIG. 1: Belhafaf et. al. experimental [^ and Monte Carlo simulation ^, |5| results on standard 
deviation of the kinetic energy distribution as a function of fragment mass number in the reaction 
235u(ni;„f). See text. 
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FIG. 2: Mass number distribution in the energy window E = 114.1 MeV for cold fission from 



233 T 



U(nift,f). Taken from Ref. [lOi]. From these data one calculates 5A= 30.14%. See text. 
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FIG. 3: Mass number distribution in the energy window 113.5 < E < 114.5 MeV for cold fission 
from ^^^U(nt/i,f). Taken from Ref. [9]. From these data one calculates 6A = 32.18%. See text. 



70 

60 

50 

g40 

a 30 

20 

10 





1 1 1 1 1 1 1 
Predicted of ■ 


- 


233,,, „ this work 


; 


♦ 


- 


- U. Quade et. al, 1988 


- 


♦ 
♦ ♦ 


- 


I.I.I. 


■ 



'80 90 100 110 

E (MeV) 



120 



FIG. 4: Odd-even effects on the proton distribution as a function of E in \]{nthS)- Values are 



taken from Ref. 



ll|, except the higher one, which corresponds to prediction of this work. See text. 
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FIG. 5: Odd-even effects on tlie neutron distribution as a function of E in U(ni/i,f). Values are 



taken from Ref. 



Ill ], except the higher one, which corresponds to prediction of this work. See text. 
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FIG. 6: Mass number distribution in the energy window 113.5 < E < 114.5 MeV for cold fission 
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from ^^U(nt;j,f). Taken from Ref. [^]. From these data one calculates 6A= 15.3%. See text. 
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FIG. 7: Mass number distribution in t 
reaction ^^^U(nt/i,f). Taken from Ref. 



le energy window E = 114.5 MeV for cold fission from the 
12l |. From these data one calculates 6A= 23.82%. See text. 
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FIG. 8: Odd-even effects on the proton distribution as a function of E in ^^^U(mf,,f). Values are 
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taken from Ref. [l^j and [13], except the higher one, which corresponds to prediction of this work. 
See text. 
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FIG. 9: Odd-even effects on tlie neutron distribution as a function of E in ^'^^U(nt/i,f). Values are 
taken from Ref. [12], except the higher one, which corresponds to prediction of this work. See text. 
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FIG. 10: Mass number distribution in the energy window E = 112 MeV for cold fission from 
^^^Pu(nt/i,f). Taken from Ref. [l5|. Prom these data one calculates 6A= 3.7%. See text. 
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FIG. 11: Mass number distribution in the energy window ii^ > 119 MeV for cold fission from 
^^^Pu(ni/i,f). Taken from Ref. [9|. From these data one calculates 5A= 7.6%. See text. 
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FIG. 12: Odd-even effects on the proton distribution as a function of E in ^^^Pu(n(ft,f). Values are 
taken from Ref. [l5[, except the higher one, which corresponds to prediction of this work. See text. 
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FIG. 13: Odd-even effects on the neutron distribution as a function of E in ^^^Pu(n(/j,f). Values 
are taken from Ref. [15[], except the higher one, which corresponds to prediction of this work. See 
text. 
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FIG. 14: Electrostatic interaction energy between spherical fragments with spherical surfaces sep- 
arated by 2 fm(C); smoothed curve of maximal available energy (Q); and maximal value of the 
total kinetic energy (K) for ^^'^U(nt/i,f). See text. 
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FIG. 15: Smoothed curves of electrostatic energy excess SC (= C — Q) and minimal values of total 
excitation energy X {= Q — K), a function of light fragment mass number A in ^^^U(n(/j,f). See 
text. 
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FIG. 16: Smoothed curves of electrostatic energy excess SC (= C — Q) and minimal values of total 
excitation energy X {= Q — K), a function of light fragment mass number A in ^'^^U(n(/j,f). See 
text. 
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FIG. 17: Smoothed curves of electrostatic energy excessiJC {= C — Q) and minimal values of total 
excitation energy X {= Q — K), a function of light fragment mass number A in ^^^Pu(ni/j,f). See 
text. 
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